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Abstract 

We describe a search for hadronic decays of the T(1S), T(2S), and T(3S) resonances to the ex- 
clusive final states p7r, K*{892)K, pa2(1320), a;/2(1270), (/>/^(1525), K*{892)K^{U30), 6i(1235)7r, 
Ki{1270)K, and Ki{UOO)K. Upper limits at 90% CL are set for all these decays from all 
three resonances below 33 X lO'^; in particular, ,B(T(1S)^ pvr) < 4 x 10"^ is the small- 
est such upper limit. For two modes, a branching fraction of zero can be ruled out with 
a statistical significance of more than 5a: B{T{1S) 0/^(1525)) = (7l^ ± 1) x IQ-^ and 
5(T(15) ^ Ki{UOO)K) = {Utt ± 2) x lO'^. Production of Ki(1270)^ in T(1S) decay is sup- 
pressed relative to that of Ki{1400)K . These results add another piece to the challenging "p — vr 
puzzle" of the charmonium system, placing constraints on models of how quantum chromodynamics 
should be applied to heavy quarkonia. All results are preliminary. 



* Submitted to the XXI International Symposium on Lepton and Photon Interactions at High Energies, 
August 2003, Fermilab, Batavia, IL, USA 
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I. INTRODUCTION AND MOTIVATION 



Little experimental information exists on exclusive decays of the T resonances below BB 
threshold. Upper limits have been published for the decays T(15') pn and T(IS') hh, 
h = n, K, p, all of the order 10-^ and T{1S) 71+71-71^ (1.8 x 10'^) ^1]. No exclusive 
final states for the T(25') and T{3S) have been examined. The situation is different in 
charmonium, where numerous channels have been measured. This by itself poses a moti- 
vation to study hadronic T decays, given the similarity of these two strongly-bound heavy 
quark systems. Furthermore, a long-standing unsolved puzzle in charmonium regarding the 
ratio of branching fractions of the 2^5*1 and l^Si (J/ip) to final states consisting of a 
pseudoscalar and a vector warrants the corresponding measurement in bottomonium. 

The expectation that the dilepton and hadronic ratios of decay widths should be at 
least roughly equal follows from QED and QCD. Both processes are thought to occur via 
annihilation of the constituent quark and antiquark, in one case to a photon and in the other 
to three gluons, and therefore are both proportional to the square of the quark- antiquark 
wave function overlap at the origin. Restating the QCD expectation for hadronic decays in 
terms of ratios of branching fractions instead of decay widths, and neglecting the running of 
the strong coupling constant^, one obtains the following prediction regarding an arbitrary 
final state H: 

_ BmS) ^ H) ^ B{^{2S) e+e-) 
^ B{J/ij ^H) ^ BiJ/ij ^ e+e-) ^ ^ 

Using the leptonic branchingfractions B{J/ip e+e") = (5.93±0.1) x 10^^ and B^ipC^S) 
e+e-) = (7.3 ± 0.4) x IQ-^ the expected value^ for the ratio is (12.3 ± 0.9)%. 

A number of channels have been studied, most of which satisfy the prediction within 
experimental errors. The most significant deviation known so far comes from the following 
vector-pseudoscalar (VP) and vector-tensor (VT) final states jl[: fm {Q < 0.007), K*{S'd2)K 
[Q < 0.001), pa2(1320) [Q < 0.02), K*{892)KI{U30) [Q < 0.02), and cj/2(1270) [Q < 
0.03 |l|,l3). 

Many theoretical approaches have been made to resolve this puzzle. None is able to 
accommodate all the measurements reported so far. For some, the crucial question is whether 
the J/ip is enhanced or the ip' is suppressed. Therefore, in addition to adding experimental 
information to the scarce amount that is available at this point on J^*^ = 1 bottomonium 
decays, in particular in the pvr channel, it is an interesting question what the ratio analog 
to Equation [T] should turn out to be for the T system. 

Depending on which theoretical model is used to explain the behavior measured in char- 
monium, the expectation for bottomonium varies. A common assumption is that the indi- 
vidual branching fractions in bottomonium should be at least one order of magnitude smaller 
than in charmonium. How much they are smaller should be of some utility in deciphering 
the puzzle. 



^ The strong coupling constant enters with the third power. The relevant ratio between the tp(2S) and 

the J/^p is {as{m^^2S)) / <^simj/^))^ = 0.85 Q- 
^ Using earlier measurements and setting the ratio of coupling constant values to unity, the ratio of dilepton 

branching resulted in a ratio Q ~ 14%. This is the reason why the puzzle posed by the failure of some 

channels was referred to as "the 14% puzzle" . 
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Two important features distinguish the situation in bottomonium from that in charmo- 
nium. First, in addition to the comparison of the 2S excitation with the ground state, 
the 3S resonance can be used due to the fact that it is below open flavor production thresh- 
old, in contrast to the situation in charmonium. Furthermore, the ratio predicted based on 
EquationQis 48% for T{2S) : T{1S) and 72% for T(3S) : T(l^). 

CLEO recently accumulated several million bottomonium decays at each of the T(IS'), 
T{2S), and T(3S') resonances. These datasets can probe decays of the J^*^ = 1 bb bound 
states at the 10~^ level. The decays pursued in this work are pvr, K*{892)K, pa2(1320), 
cu/2(1270), 0/^(1525), K*{892)K*{U30), 6i(1235)7r, Ki{1270)K, and Ki{UOO)K. They 
sample VP, axialvector-pseudoscalar (AP), and VT type final states with and without 
strangeness and constitute the most copious two-body hadronic final states in J/ip decay 
(each with a branching ratio of 0.1 — 1%). Each proceeds via the strong interaction and 
conserves isospin. These decay modes provide a logical starting point for using the T system 
to untangle the many subleties of the pvr puzzle and associated anomalies. 

II. ANALYSIS OVERVIEW 

The analysis strategy is straightforward. Event selection criteria for the different modes 
are developed using signal Monte Carlo. We emphasize cleanliness over efficiency in order 
to suppress decays faking the desired final states. CLEO data taken at and just below the 
T(4S') resonance, suitably scaled by luminosity and beam energy, is used as an indication of 
background properties and final contamination levels. Finally, projected backgrounds and 
event totals are normalized by efficiencies and the number of resonance decays produced, 
and branching fractions or upper limits computed. 

It is important to note that many different kinds of backgrounds contribute to the sample 
obtained at the T(4S') resonance. Not all of them scale with the same center-of-mass energy 
dependence (see discussion below). For that portion of backgrounds which are truly the 
same final state but produced non-resonantly, that is, proceed as e+e~ — 7* — > if instead 
of e'^e^ ^ T H, there is also the possibility of interference, which has been neglected in 
this work. 

III. EVENT SELECTION 

The CLEO III detector is described in detail in IQ]. Its key features exploited in this 
analysis are a solid angle coverage for charged and neutral particles of 93% and two par- 
ticle identification systems to separate kaons from pioiis, namely using energy loss in the 
drift chamber and a Ring Imaging Cherenkov detector [7|. The tracking system achieves a 
charged particle momentum resolution of 0.35% (1%) at p =1 GeV/c {p =5 GeV/c) and the 
calorimeter a photon energy resolution of 2.2% (1.5%) at E-y =1 GeV {E^ =5 GeV). The 
combined dE/dx-KlCB. particle identification system attains a kaon efficiency (fake rate) 
>90% (< 5%) below 2.5 GeV/c and falls (rises) to -70% (-25%) near 5 GeV/c. 

Standard requirements are used to identify charged particles from tracks in the drift 
chamber and photons from electromagnetic showers in the Csl calorimeter. The simplicity 
of the final states under study is exploited by imposing an energy conservation requirement 
on Xt = T,h{Eh) / {2Ebeam) of 0.98 < Xt < 1.015 when summing over the energies of the 
decay products. For all the target modes, the experimental resolution in this quantity is 
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smaller than 1%. For each of the final state resonances, a search window for the invariant 
mass of the decay products is established based on signal Monte Carlo studies. Electron 
and muon vetoes are imposed to suppress QED backgrounds. 



IV. MONTE CARLO SAMPLES 

Efficiencies are evaluated with Monte Carlo simulation of the process and detector re- 
sponse [sf. The pvr and K*K modes are generated with the polar angle 9 distributed accord- 
ing to (1 + cos^O). All other channels are thrown to be fiat in cos^ as they can be of any 
linear combination of sin^ 6 and (1 + cos^ 6) . The final efficiencies are of the order of 5 — 10%, 
including all effects of selection criteria and all intermediate branching fractions (Table P). 
The detection efficiency of some of the modes varies significantly with the beam energy and 
so is evaluated at the T(IS') and T(45') separately. The T{2S) and T(3S') efficiencies, also 
listed in Table lU are obtained by interpolating linearly between the T(IS') and T(4S'). 

Systematic uncertainties on the efficiency include uncertainties in the polar angle distri- 
bution for VT and AP modes (5%) and modeling of tracks (1% per charged track), vr^s (8% 
per beam-energy- tt" and 5% per softer 7r°), lepton veto (1% per particle), kaon identification 
(5% per identified kaon) and pion fake rate (3%), and secondary vertex- finding (5%). These 
contributions and the uncertainty in the number of produced T decays (5%) are summed in 
quadrature. The resulting total relative error is close to 10% for all modes. 

Although r-pair production of the states in question contributes for Xt < 0.98, the 
background in the signal region is found to be small, as is that from /i-pairs. 

V. DETERMINATION OF SIGNAL AND BACKGROUND YIELD 

The data samples used consist of (21 ±1), (5.4±0.2), and (5.0±0.3) xlO^ T(1S), T(2^), 
and T(3S') decays, respectively.^ Events that satisfy the event selection criteria mentioned in 
the previous section are simply counted. To obtain maximal statistical power, both charged 
and neutral decay modes were considered and combined for the final result. The uncorrected 
event numbers thus obtained are listed in Table UTl Example distributions of the scaled total 
energy and invariant masses for the nine channels under study for T(IS') decays can be found 
in Figures n and 121 

The background level was determined using a large amount of data taken on or near 
the T(45') resonance, extrapolated down to the lower resonances. For the extrapolation, 
three issues must be taken into account: Scaling with luminosity, efficiency dependence 
on the center-of-mass energy, and cross section dependence on the center-of-mass energy. 
The efficiency dependence can be read off Table lU The ratio of integrated luminosities 
for the datasets used is 0.133, 0.093, and 0.135 for T{IS) : T(45), T(25) : T(45), and 
T(3S') : T(4S'), respectively. The cross section extrapolation with beam energy poses the 
most uncertain contribution. Depending on the contributing background process, it could 
vary from 1/s for QED processes to as much as l/s^ (l/'S^ for channels such as pii or 



K*K).[^ We therefore quote ranges of scale factors that take this uncertainty into account. 



shown in Table ITTTl The scale factors vary between 0.10 and 0.45. Within a specific channel, 



^ The number of resonance decays was computed using preliminary hadronic cross section line shape mea- 
surements with CLEO III by scanning the center-of-mass energy around the resonance. 
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the ratio of upper to lower scale factor limit ranges from a maximum of almost 2 for the 
T(IS') VP modes down to the level of 1.1 for T{3S) channels. 

Cross-feed between the investigated final states is accounted for as a separate source of 
background. The largest contribution is found to be i^i(1400)^ events leaking into the 
Ki{1270)K sample. The Ki{UOO)K signal above M{Knn) = 1.4 GeV (where there is 
no i^i(1270) signal) is evaluated with only the T(4S') background considered and then, 
scaled in accordance to the Monte Carlo detection efficiency, taken as a second source of 
background for the Ki{1270)K. We neglect cross-feed in the other direction because the 
predicted backgrounds to Ki{1270)K saturate the observed rate. Also, the pn channel 
receives contributions from mis-identified K*{892)K events in the p~7i~^ mode. Again, the 
K*{892)K leakage is treated as a second source of background, properly scaled. Using 
the T(4S') event yield from Table ITTl together with the maximal scale factors from Table ITTTl 
and, where necessary, scaled cross-feed contributions, one arrives at estimates of background 
levels listed in Table HVl 

The confidence level that any given mean signal combined with background would exceed 
or equal the observed event count is computed from simulated trials in which a pseudo- 
random number generator is employed to throw Poisson distributions. Poisson fiuctuations 
in both the observed resonance and 4S samples are simulated by allowing not only the 
background to vary around its mean from one trial to the next, but also the mean background 
itself: 4S levels are fiuctuated around the observed number prior to application of the scale 
factor to obtain the mean, and only then are fiuctuations on the mean introduced. We reject 
thrown backgrounds which exceed the number of observed events. Trials are thrown in steps 
of 0.1 in signal mean until the desired confidence level is exceeded. This procedure predicts 
slightly larger intervals than the approach by Feldman and Cousins ^ when backgrounds 
are less than the observed number of events, and considerably wider ones for observations 
smaller than the mean expected background. 

Upper limits at 90% CL on the number of signal events are listed in Table for which 
the 1/s scale factors were used to minimize the estimated non- resonant background and 
therefore to maximize any potential signal. These are converted into 90% CL upper limits 
on the corresponding branching fractions shown in Table IVll using the number of produced 
resonance decays and efficiencies (Table EJ; we account for the systematic relative error of 
10% in this conversion by increasing each upper limit by an additional 1.28o" = 13%. 

Two-sided intervals of 68% confidence level are also shown in Table IVTl for channels with 
statistical significance exceeding one standard deviation. We define the statistical signifi- 
cance to be the number of Gaussian standard deviations above which lies the probability that 
the background alone fiuctuated up to the observed number of events. For these intervals, 
the upper end of the 4S scale factor range is employed to minimize the chance of undersub- 
tracting background. The systematic error shown includes the uncertainty on the product of 
efficiency and number of T's produced mentioned above (10%) in quadrature with an addi- 
tional 10% to account for uncertainties in e~^e~ annihilation backgrounds as well as T cross- 
feed estimates. For two channels, B{T{1S) 0/^(1525)) and B{T{IS) Ki{im)K), 
zero branching fraction can be ruled out at a statistical significance of more than 5a, mak- 
ing these the first exclusive hadronic decay modes measured in the T system. Several other 
channels with significance near ~ 3a show suggestive but statistically marginal evidence for 
branching fractions at the few per million level. 

In contrast to the ip{2S) and in similarity to the J/ip [lo!], the production of Ki{1270)K 
is suppressed relative to that of i^i(1400)^. 
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VI. CONCLUSIONS 



Using CLEO III datasets of 21, 5.4, and 5.0 million T(1S), T(2S), and T(3S) decays, 
respectively, we have searched for nine of the potentially most probable two-body-hadronic 
decays, which include VP, AP, and VT channels: jm, K*{892)K, pa2(1320), a;/2(1270), 
(/)/^(1525), ir* (892) is:* (1430), 6i(1235)7r, Ki{1270)K, and Ki{UOO)K. The upper limit 
at 90% confidence level for T(1S)— > pre is lowered by more than an order of magnitude 
to 4 X 10^^, and 90% CL upper limits for the other eight modes, measured for the first 
time, range from 7 — 33 x 10^®. Two channels have been observed at convincing statistical 
significance: S(T(1S) ^ 0/^(1525)) = {7tl ± 1) x lO'^ and B{T{1S) ^ i^i(1400)i^) = 
(14^3 ±2) X 10~^. The branching fractions from the T(1S) are smaller than the comparable 
values on the J/ip by factors of several hundred (for i^i(1400)A') to at least several thousand 
(for pvr); the T(15') — > pvr branching fraction is measured to be suppressed by at least six 
powers of Mj/^/My relative to J/ip ^ pit. The above results are preliminary. 
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FIG. 1: Scaled total energy, Xt, distribution for candidate channels from T(1S), in which all 
selection criteria except the Xt requirement have been applied. Solid points represent data; the 
red (dark) histogram is signal MC with arbitrary normalization; the green (light) histogram is 
scaled 4S data. 
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FIG. 2: Mass distributions of intermediate states Mx, with X given by the label in parentheses, 
for candidate channels from T(1S) decays, in which all selection criteria except the one on the 
plotted mass have been applied. Solid points represent data; the red (dark) histogram is signal 
MC with arbitrary normalization; the green (light) histogram is scaled 4S data. Plots on the right 
refer to the same channel as each neighboring plot on the left. 
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TABLE I: Efficiency in percent for the isospin- and charge-conjugate- inclusive parent decay chan- 
nels listed, summed over the various sub-modes used in this analysis, including all the effects of 
selection criteria and all intermediate branching fractions. 



QiMUuel 


T(1S) 


T(2S) 


T(3S) 


-IS 




7.8 


6.6 


6.0 


5.6 


K*{892) K 


10.6 


10.0 


9.6 


9.5 


p 02(1320) 


7.8 


7.0 


6.6 


6.3 


u /2(1270) 


7.7 


6.9 


6.5 


6.1 


^ /^(1535) 


10.1 


10.1 


10.2 


10.2 


K*{m2) ii'|(1430) 


5.3 


5.2 


5.1 


5.1 


61(1235) TT 


7.2 


6.6 


6.2 


6.0 


is:i(1270) K 


9.2 


8.9 


8.7 


8.6 


i^i(1400) K 


9.4 


9.3 


9.3 


9.3 



TABLE IL Number of events in each isospin- and charge-conjugate-inclusive decay channel. 



Channel 


T(1S) 


T(2S) 


T(3S) 


4S 


P TT 


4 


1 


3 


6 


A'* (892) A' 


18 


2 


1 


15 


p 02(1320) 


29 


8 


10 


47 


CO /2(1270) 


6 


1 





4 


cP /^(1535) 


17 


4 





7 


a:* (892) A'|(1430) 


16 


6 


5 


23 


61 (1235) TT 


6 


1 


2 


3 


A:i(1270) K 


27 


7 


8 


29 


A'i(1400) K 


37 


13 


9 


38 
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TABLE III: Range of scale factors for the 4S yield to the lower resonances for the combined isospin- 
and charge-conjugate-inclusive decay channels, including luminosity, efficiency energy dependence, 
and range of l/s" scaling (see text). 



Channel 




Tf2S) 


T(3S) 


p TT 


0.23-0.45 


0.12-0.17 


0.15-0.17 


K*(892) K 


0.19-0.36 


0.13-0.15 


0.14-0.16 


p 02(1320) 


0.21-0.32 


0.11-0.14 


0.15-0.16 


CO /2(1270) 


0.21-0.32 


0.12-0.14 


0.15-0.16 


</> /2(1535) 


0.17-0.26 


0.12-0.13 


0.14-0.15 


K*{892) 7^:1(1430) 


0.17-0.27 


0.12-0.13 


0.14-0.15 


6i(1235) TT 


0.20-0.31 


0.11-0.14 


0.15-0.16 


Ki(1270) K 


0.17-0.27 


0.12-0.13 


0.14-0.15 


Ki(1400) K 


0.16-0.25 


0.12-0.13 


0.14-0.15 



TABLE IV: Approximate background level estimate(s) from non-resonant e~^e~^— annihilation 
scaled from the 4S data, and, where indicated by a "+", cross-feed from other channels (see text), 
in the isospin- and charge-conjugate- inclusive decay channels. 



Channel 


T(1S) 


T(2S) 


T(3S) 


P TT 


3+2 


1+0 


1+0 


K*{892) K 


5 


2 


2 


p 02(1320) 


15 


7 


8 


CO /2(1270) 


1 


0.6 


0.6 


/2(1525) 


2 


1 


1 


i^*(892) K2*(1430) 


6 


3 


4 


fei(1235) TT 


1 


0.4 


0.5 


ii:i(1270) K 


8+20 


4+8 


4+4 


^1(1400) K 


10 


5 


6 



11 



TABLE V: Upper limits on number of signal events seen in each channel, at 90% confidence level, 
computed using the observed 4S yield with the lower end of the range of scale factors (i.e. based 
on a 1/s cross-section dependence, resulting in a low background prediction), and including the 
statistical fluctuations of not only signal and background but also the mean of the background. 



Channel 


T(1S) 


T(2S) 


T(3S) 


p TT 


5.5 


3.5 


5.9 


if* (892) K 


22.1 


4.0 


6.1 


p 02(1320) 


27.6 


8.1 


8.8 


/2(1270) 


9.8 


3.6 


2.3 


</> /2(1535) 


22.5 


7.3 


2.3 


X*(892) K|(1430) 


18.7 


8.0 


6.4 


6i(1235) TT 


10.0 


3.7 


5.0 


Ki{m{)) K 


13.3 


4.6 


6.7 


ii:i(1400) K 


40 


14.5 


9.2 



TABLE VI: T branching fraction 68% confidence intervals and 90% CL upper limits, in units of 
10~^, for each isospin- and charge-conjugate-inclusive decay channel. The first error is statistical 
and the second is systematic. For the 90% CL upper limits, a low-end 4S scale factor was used, 
whereas for the 68% confidence intervals the maximal scale factor value was employed. Also shown 
is the statistical signal significance (in units of a). A "— " indicates less than la significance, 
defined as the number of Gaussian standard deviations above which lies the probability that the 
background alone fluctuated up to the observed number of events. 



Channel 


T(1S) 




T(2S) 




T(3S) 






Interval/Sig. 


UL 


Interval/Sig. 


UL 


Interval/Sig. 


UL 


P TT 




4 




11 


9i| ± 1 / 1.5 


22 


if* (892) K 


6l^ ± 1 / 3.6 


11 




8 




14 


p a2(1320) 


9±4 ± 1 / 3.0 


19 




24 


±1/1.3 


30 


a; /2(1270) 


3ii ± 1 / 2.6 


7 




11 




8 


4> /2(1525) 


7^^ ± 1 / 5.5 


12 


6^3 ± 1 / 3.0 


17 




14 


K*(892) K^{Wi{)) 


9^:^ ± 1 / 3.0 


19 


11±8±2 / 1.6 


32 




28 


6i(1235) TT 


3±2 ±1/2.9 


8 




12 


5l| ± 1 / 1.4 


18 


ii:i(1270) K 




8 




11 




17 


Ki(1400) K 


±2 / 5.(3 


23 


IGl"^? ± 2 / 2.9 


33 


1^1 ±1/1.5 


22 



12 



